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BLENDING  METHODOLOGY  OF  LINEAR  PARAMETER  VARYING  CONTROL 
SYNTHESIS  OF  F-16  AIRCRAFT  SYSTEM 


JONG-YEOB  SHIN*,  GARY  J.  BALAS+,  AND  ALPAY  M.  KAYAi 

Abstract.  This  paper  presents  the  design  of  a  linear  parameter  varying  (LPV)  controller  for  the  F-16 
longitudinal  axes  over  the  entire  flight  envelope  using  a  blending  methodology  which  lets  an  LPV  controller 
preserve  performance  level  over  each  parameter  subspace  and  reduces  computational  costs  for  synthesizing 
an  LPV  controller.  Three  blending  LPV  controller  synthesis  methodologies  are  applied  to  control  F-16 
longitudinal  axes.  Using  a  function  substitution  method,  a  quasi-LPV  model  of  the  F-16  longitudinal  axes 
is  constructed  from  the  nonlinear  equations  of  motion  over  the  entire  flight  envelope,  including  non-trim 
regions,  to  facilitate  synthesis  of  LPV  controllers  for  the  F-16  aircraft.  The  nonlinear  simulations  of  the 
blended  LPV  controller  show  that  the  desired  performance  and  robustness  objectives  are  achieved  across  all 
altitude  variations. 

Key  words.  LPV  control  synthesis,  F-16  longitudinal  axes 
Subject  classification.  Guidance  and  Control 
Nomenclature. 

TO,  lyy  :  Mass  (slug).  Inertial  moments  (slug  ft^) 

V,q  :  Velocity  (ft  sec“^).  Dynamic  pressure  (psi) 

a,  q  :  Angle  of  attack  (rad).  Pitch  rate  (rad  sec“^) 

0,  7  :  Pitch  angle  (rad).  Flight  path  angle  (rad) 

5e,T  :  Elevator  deflection  (rad).  Thrust  (lb) 

c,  S  :  Chord  length  (ft).  Reference  area  (ft^) 

Xac,  Xcg  ■  Aerodynamic  center  position  (ft).  Center  of  gravity  position  (ft) 

Cx,  Cz  ■  X  and  Z  force  aerodynamic  coefficients 

Cxg ,  Czg  ■  Aerodynamic  stability  derivatives 

Cm„  ■  Pitch  moment  aerodynamic  coefficient 

Crriq  :  Pitch  moment  aerodynamic  stability  derivative 

1.  Introduction.  Extensive  researches  have  focused  over  the  last  ten  years  on  developing  analysis 
and  synthesis  techniques  for  gain-scheduled  controllers  for  linear  parameter  varying  (LPV)  systems  [13,  12, 
17,  6,  2,  15].  In  Ref.  [13,  12],  conditions  are  given  which  guarantee  stability,  robustness,  and  performance 
properties  of  the  global  gain-scheduled  designs.  Recent  theoretical  developments  have  produced  methods 
for  synthesizing  gain-scheduling  controllers  for  LPV  systems,  which  guarantee  a  level  of  robust  stability 
and  performance  across  scheduling  parameter  spaces  [6,  2,  15].  In  Ref.  [2],  LPV  control  synthesis  methods 
have  also  been  developed  using  parameter-dependent  Lyapunov  functions  to  lead  to  a  less  conservative 
result.  This  gain-scheduling  approach  has  been  successfully  applied  to  synthesize  controllers  for  pitch- axis 
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missile  autopilots  [3,  11],  F-14  aircraft  lateral-directional  axis  during  powered  approach  [4,  5],  and  turbofan 
engines  [8]. 

One  of  the  potential  difficulties  in  practical  uses  of  the  LPV  synthesis  methodology  with  parameter- 
dependent  Lyapunov  functions  is  that  the  complexity  of  the  linear  matrix  inequality  (LMI)  optimization 
problem  increases  exponentially  with  the  number  of  scheduling  parameters  and  the  number  of  the  grid 
points  over  the  scheduling  parameter  spaces.  One  approach  to  facilitate  practical  use  of  LPV  synthesis 
methodology,  denoted  as  the  “blending  approach”,  has  been  discussed  in  Ref.  [15,  14].  This  approach  to 
control  design  partitions  the  full  parameter  space  into  overlapping  small  subspaces.  An  LPV  controller  is 
synthesized  for  each  small  region.  These  regional  controllers  are  blended  into  a  single  LPV  controller  for  the 
entire  parameter  space.  In  this  paper,  this  blending  approach  is  applied  to  control  F-16  longitudinal  axes 
over  the  entire  flight  envelope. 

To  synthesize  an  LPV  controller  for  an  F-16  aircraft,  an  LPV  model  of  the  aircraft  dynamics  is  required. 
Conventional  approaches  to  generate  an  LPV  model  of  an  aircraft  are  based  on  Jacobian  linearizations  at 
trim  points  or  a  change  of  state  coordinates  [11]  to  reduce  the  nonlinearity  of  aircraft  dynamics.  The  LPV 
models  constructed  by  both  approaches  can  present  aircraft  dynamics  at  trim  conditions.  However,  the 
models  can  not  represent  aircraft  dynamics  at  non-trim  conditions.  Instead  of  using  Jacobian  linearization 
or  state  transformation,  the  nonlinear  terms  of  aircraft  dynamics  can  be  substituted  for  other  functions  in 
quasi-LPV  form  [18,  19].  This  function  substitution  approach  can  be  applied  for  both  trim  and  non-trim 
conditions.  The  approach  has  been  used  in  generating  a  quasi-LPV  model  of  a  generic  missile  [18,  19].  In 
this  paper,  a  quasi-LPV  model  of  F-16  longitudinal  axes  is  provided  over  the  entire  flight  envelope  including 
non-trim  regions,  using  a  function  substitution  approach. 

In  Section  2,  a  brief  summary  of  conventional  LPV  controller  synthesis  used  in  this  paper  is  presented 
to  emphasize  the  complexity  of  the  LMI  optimizations.  In  Section  3,  three  blending  LPV  control  synthesis 
methodologies  are  presented.  Development  of  a  quasi-LPV  model  of  F-16  longitudinal  axes  is  presented  in 
Section  4.  In  Section  5,  formulation  of  the  LPV  control  problem  and  the  blending  of  two  LPV  controllers  of 
the  F-16  aircraft  are  presented.  Nonlinear  simulations  of  the  closed-loop  system  with  the  blended  controller 
are  presented  in  Section  6  and  this  paper  concludes  with  a  brief  summary  in  Section  7. 

2.  LPV  Control  Synthesis.  In  this  section,  a  conventional  LPV  control  synthesis  using  parameter- 
dependent  Lyapunov  functions[2]  is  briefly  described.  Consider  a  generalized  linear  open-loop  system  as 
functions  of  parameters  p{t)  €  V.  For  a  compact  subset  V  C  TZ^ ,  the  parameter  variation  set  denotes  the 
set  of  all  piecewise  continuous  functions  mapping  R  (time)  into  V  with  a  finite  number  of  discontinuities  in 
any  interval,  where  s  is  the  number  of  parameters.  An  LPV  open-loop  system  can  be  written  as 

x{t)  Mp{t))  Bi{p{t))  B2{p{t))  x{t) 

e{t)  =  Ci{p{t))  0  Di2(p(t))  d{t)  ,  (2.1) 

C2{p{t))  D2i{p{t))  0  _  u{t)_ 

where  y{t),  e{t),  d{t)  and  u{t)  are  measurements,  errors,  disturbances,  and  control  signals.  Hereafter,  p 
denotes  p{t).  The  induced  C2  norm  of  d  to  e  is  defined  as 

l|e||2 

TUFT- 

p€V,d€C-i,\\d\\ii=0  ll“l|2 

In  a  conventional  LPV  synthesis  methodology,  suppose  there  is  an  LPV  output  feedback  controller  K  (p) 
which  stabilizes  the  closed-loop  system  exponentially  and  makes  the  induced  T2-norm  of  d  to  e  less  than  7. 
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The  controller  K  (p)  can  be  written  as 


Xk 

Ak{p) 

Bkip) 

Xk 

u 

Ck{p) 

Dk{p) 

y 

(2.2) 


An  LP V  controller  K  (p)  can  be  constructed  from  solutions  of  X  (p)  and  Y  (p)  of  the  following  optimization 
problem  [2]: 


subject  to 


mm  7 
x,Yen”'^” 


^  a  Y 

X{p)AT{p)+A{p)X{p)-Y,{^i^)-B2{p)Bl{p)  X{p)CUp)  7-'Si(p) 

Cii{p)X{p) 

l-^Bfip)  0 


0 

-In 


<0, 


(2.3) 


(2.4) 


A{p)Y{p)  +  Y{p)AT{p)  +  ^(£,^)  -  Cj{p)CAp)  Y{p)B,Ap)  l-^Cjip) 

i=l 


bUp)y{p) 

l~^Ci{p) 


-I, 


ndi 

0 


0 

-In 


<  0, 


(2.5) 


X(p)  7-'4  ,  „ 

_7-'4  Y{p)  \  -  ’ 

X(p)>0,  Y(p)>0, 

where 

A{p)  =  A{p)  -  BAp)C,Ap),  Mp)  =  A{p)  -  B,Ap)C2{p), 


(2.6) 


(2.7) 


and  n  is  the  number  of  states  of  the  generalized  open-loop  system.  The  detailed  definitions  of  matrices 
C'ii(p),  C'i2(p),  and  Bii(p)  can  be  found  in  Ref.  [2].  Note  that  indicates  that  every  combination  of 

i=l 

and  Pi  is  included  in  the  LMIs.  The  parameter  rate  p  is  bounded  as  Hi  <  Pi  <  Pi- 

A  method  to  construct  an  LPV  controller  K (p)  from  the  solutions  X  (p)  and  Y  (p)  is  taken  from  Ref.  [2] . 


In  this  paper,  an  LPV  controller  is  constructed  as  [2]: 

Ak{p)  =  A{p)  +  B2{p)F{p)  +Q-\p)Y{p)L{p)C2{p)  -7-^Q-\p)M{p,p),  (2.8) 

Bk{p)  =  -Q-Hp)Y{p)L{p),  (2.9) 

Ck{p)=F{p),  (2.10) 

Dk(p)=0,  (2.11) 

where  matrices  Q{p),  F{p),  L{p),  and  M{p,p)  are  defined  as 
Q(p)=Y(p)-j-^X-^(p), 

F{p)  =  -[Bl{p)X-Hp)  +  DUp)Ci{p)], 

L{p)  =  -[Y-Ap)c^{p)  +  bAp)dUp)], 


M{p,p)  =  H{p,p)  +^‘^Q{p)[-Q-^{p)Y{p)L{p)DM  -  Bi{p)]B^{p)X-^{p). 
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Matrix  H{p,p)  is  defined  as 


H{p,p) 


^  \ 

-  [X-\p)Ap{p)+Ap{pfX-\p)+X[Pi^) 

+  Cl  {p)Cf  {p)  +  7-2X-1  (p)Bi  (p)Bi  {pfX-^  (p)] 


with  Ap{p)  =  A{p)+B2{p)F{p)  and  Cf{p)  =  Ci{p)+Di2{p)F{p).  The  closed-loop  system  with  the  controller 
K{p)  is  exponentially  stable  and  the  induced  C2  norm  is  less  than  7.  The  proof  can  be  found  in  Ref.  [2]. 

To  make  the  optimization  problem  of  equation  (2.3)  computationally  tractable,  the  scheduling  param¬ 
eters  p  are  discretized  into  grid  points.  Thus,  an  infinite  number  of  constraints  are  represented  by  a  finite 
number  of  LMI  constraints.  Also,  X{p)  and  Y (p)  are  represented  by  a  finite  number  of  basis  functions  hi{p) 
and  gi{p): 

Ny 

x{p)  =  E 

where  hi{p)  and  gi{p)  are  continuously  differentiable  functions.  The  LMIs  of  equations  (2.4)-(2.6)  are  solved 
for  all  grid  points  of  the  scheduling  parameters  simultaneously.  The  size  of  the  optimization  problem  is 
proportional  to  2^~^^Ng  where  s  and  Ng  are  the  number  of  scheduling  parameters  and  total  number  of  grid 
points  over  the  scheduling  parameter  space.  The  number  of  decision  variables  are  +  Ny)  ”(”+^)  ^  where 
Nx,  Ny  and  n  are  the  order  of  basis  functions  of  X  and  Y,  and  the  state  order.  Thus,  computational  time 
to  solve  the  optimization  problem  of  equation  (2.3)  is  dependent  on  the  number  of  grid  points  of  scheduling 
parameters,  the  state  order  of  a  generalized  open-loop  system,  and  basis  function  orders.  The  conventional 
LPV  synthesis  methodology  may  require  expensive  computational  cost  (computational  time  and  computer 
memory)  when  number  of  scheduling  parameters  increases.  In  the  next  section,  a  blending  method  to  reduce 
the  computational  cost  to  synthesize  an  LPV  controller  and  to  preserve  performance  level  over  parameter 
subspaces  is  presented. 

3.  LPV  Controller  Blending  Approach.  Instead  of  designing  a  single  LPV  controller  for  the  entire 
parameter  set  in  a  conventional  LPV  synthesis,  LPV  controllers  can  be  synthesized  for  parameter  subsets, 
which  are  overlapped  with  each  other.  Then,  these  LPV  controllers  are  blended  into  a  single  LPV  controller 
over  the  entire  parameter  set.  Thus,  the  performance  of  the  closed-loop  system  with  the  blended  controller 
is  preserved  when  parameter  trajectories  travel  over  the  overlapped  parameter  subsets. 

Consider  the  scheduling  parameter  vector  p  €  P  which  consists  of  subvectors  ps  €  Vs  and  Pd  &Vd-  The 
parameter  subset  Vs  can  be  partitioned  into  two  subsets  with  the  following  conditions: 


Vs^=Vs,r)Vs„  (3.1) 

Vs  =  Vs,  U  Vs,.  (3.2) 


Suppose  that  there  exist  LPV  controllers  Ki  and  K2  constructed  from  parameter-dependent  Lyapunov 
functions  of  a  parameter  subvector  pd,  over  each  parameter  subset  Vs,  U  Vd,  i  =  'i-,  2.  Thus,  V,  and  1) 
are  functions  of  pd  (not  ps)  over  each  parameter  subset  Vs,  U  VdA  =  1,  2.  Also,  the  controller  A,  can 
stabilize  the  closed-loop  system  and  the  induced-T2  norm  of  the  closed-loop  system  is  less  than  7,  over  each 
parameter  subset  Vs,  U  Vd- 

When  a  scheduling  parameter  subvector  ps  is  in  the  intersection  subset  Vs,, ,  LMI  solution  matrices 
Xi{pd)  and  Yi(pd)  are  combined  into  Xh{ps,pd)  and  Yh{ps,pd),  respectively.  There  are  three  methods  to 
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calculate  Xf,lyPs,pd)  and  Yf,lyPs,pd),  which  will  be  explained  later.  The  blended  matrices  Xf,lyPs,pd)  and 
Yb{ps,Pd)  should  be  feasible  solutions  of  the  LMI  constraints  of  equations  (2.4)-(2.5)  over  the  parameter 
subset  Vs^  U  Vd-  When  the  parameter  subvector  ps  is  in  the  parameter  subset  Vsi  —  Xb{ps,Pd)  and 
Yb{ps,Pd)  should  be  equal  to  Xi{pd)  and  Yi{pd),  respectively,  and  the  partial  derivatives  of  and 
dX  ■  dY- 

should  be  equal  to  and  respectively.  There  are  three  blending  methods  for  Xb  and  Yb  to  satisfy 
the  feasibility  condition  and  the  boundary  conditions. 

For  method  I,  matrices  Xb(ps,Pd)  and  Yb(ps,pd)  can  be  written  as: 

2 

XbiPs,  Pd)  =  hxi  {ps)Xi{pd),  (3.3) 


i=l 


Yb{ps,Pd)  =  '^KiPs)YiiPd), 


(3.4) 


i=l 


where  “blending  functions”  bx^ips)  and  by.{ps)  are  differentiable  scalar  functions.  To  satisfy  the  boundary 
conditions  of  Xb  and  Yb,  blending  functions  bx^ips)  and  by.{ps)  are  defined  as: 


bxi (Ps)  — 


1,  Ps  G  Vsi  —  Vsr,  , 

0,  Ps  €  Vs  —  Vsi  , 


(3.5) 


iPs)  — 


1,  Ps  G  Vsi  —  Vs 
0,  Ps  €  Vs  —  Vsi 


dbx(Ps)  „  dby.(ps)  „  ^  ^ 

’  =  0,  =0,  psG  Vs,  -  Vs^ . 

dps  dps 

Suppose  that  the  blending  functions  bxi(ps)  and  by.(ps)  satisfy  the  following  conditions: 

2 

'^bxi(ps)  = 0  <bxi(ps)  ps^Vsr,, 

2 


by,  (Ps)  =  1,  0  <by,  (Ps)  <1,  Ps  &Vs 


(3.6) 

(3.7) 

(3.8) 

(3.9) 


i=l 


and  7  is  chosen  as  max(7i,72).  Then,  the  following  equations  are  satisfied,  since  the  LMIs  of  equations 
(2.4)-(2.5)  are  convex  with  respect  to  X  and  Y.  Hereafter,  p  dependence  is  omitted  for  convenience. 

^  ^  dX 

XbA^  +  AXb  -  W(£,  W  bx,  ^)  -  B2B[  +  XbC[iCiiXb  +  ^-^B,Bf  <  0,  (3.10) 

^  yt 

^  ^  dY- 

AYb  +  YbA^  +  Y(£i  E  +  YbBuB^^Yb  +  <  0.  (3.11) 

i=i  j=i  P* 

db  (  1  db  (  1 

When  the  derivatives  of  blending  functions  — gp  and  — Qp  are  small  enough  to  satisfy  the  following 

inequalities. 


1  <mx), 

[  i=i 


dpi 


7=1  I  i=i 


dpi 


(3.12) 

(3.13) 
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it  is  obviously  noted  that  the  blended  Xi,{ps,pd)  and  Yh{ps,pd)  can  be  feasible  solutions  of  LMIs  of  equations 
(2.4)-(2.6)  over  the  parameter  subset  Vs^ ,  where  matrices  Mx  and  My  are  defined  as 


dXi 


Mx  =XkA^  +  AX,  - 


s  2 


dYy 


My  ^AY,  +  Y,A^  +  ^(£,;  ^  hy, 

i=l  7=1 


B^B'^  +  X,Cf^CuX,  +  , 


c^c2  +  nsiiBfin  + 


Here,  a  and  a  represent  the  maximum  and  the  minimum  singular  values,  respectively. 

The  procedure  of  designing  an  LPV  controller  over  the  entire  parameter  set  using  this  blending  method 

I  is: 

1.  Partition  the  entire  parameter  set  into  two  subsets  which  have  the  overlapped  subset. 

2.  Solve  the  LMI  optimization  of  equation  (2.3)  for  Xi{pd)  and  Y^pd)  over  each  parameter  subset. 

3.  Define  blending  scalar  functions  hxAps)  and  hyAps)  which  satisfy  the  boundary  conditions  of  equa¬ 
tions  (3.5)-(3.7)  and  the  derivative  conditions  of  equations  (3.12)-(3.13).  Note  that  the  derivative 
conditions  are  sufficient  for  X,{ps,pd)  and  Y,{ps,pd)  to  be  feasible  solutions  of  the  LMI  constraints 
of  equations  (2.4)-(2.5)  over  the  parameter  subset  Vst^  A  controller  designer  chooses  candidate 
blending  functions  until  the  feasibility  conditions  of  equations  (2.4)-(2.5)  are  satisfied. 

4.  Construct  an  LPV  controller  over  the  entire  parameter  set  VsUVd  using  equations  (2.8)-(2.11),  from 
the  calculated  X,  and  Y,  of  equations  (3.3)-(3.4). 

For  method  II,  an  alternative  way  to  calculate  the  blended  matrices  X{ps,pd)  and  Y{ps,pd)  as  feasible 
solutions  of  the  LMIs  of  equations  (2.4)-(2.6)  is  adding  extra  LMIs  in  conventional  LPV  synthesis  with 
candidate  blending  functions  [2].  The  extra  LMIs  are: 

-  2  2  s  „  2 

^  bx.x^A^  +  i ^  bx,x^  +  E 

i=l  i=l  i=l  i=l 

2 

i=i 

7-1 

^  E  K  Yj  +  E  YjA^  +  E(£.  ^  E  K  Yj)  -  cl 

i=l  i=l  i=l  i=l 

2 

i=i 

7“'C^i 


Y^byXjBu  1-^Cj 

i=i 


0 


<  0. 


(3.15) 


Y^hx.XfiA  7-IB1 

i=i 


-4, 

0 


-4 


<  0, 


(3.14) 


The  procedure  of  designing  an  LPV  controller  over  the  entire  parameter  set  using  this  blending  method 
II  is: 

1.  Partition  the  entire  parameter  set  into  two  subsets  which  have  the  overlapped  subset. 

2.  Solve  the  LMI  optimization  of  equation  (2.3)  for  the  solution  matrices  Xi{pd)  and  Yi{pd)  over  one 
of  parameter  subsets,  which  is  denoted  by  Vs^  U  Vd- 

3.  Define  blending  scalar  functions  bxAps)  and  byAps)  which  satisfy  the  boundary  conditions  of  equa¬ 
tions  (3.5)-(3.7). 
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4.  Solve  the  LMI  optimization  problem  of  equation  (2.3)  with  the  extra  LMI  constraints  of  equa¬ 
tions  (3.14)-(3.15)  for  the  solution  matrices  X^ipd)  and  Y^ipd)  over  the  other  parameter  subset 
Vs,  U  Vd- 

5.  Construct  an  LPV  controller  over  the  entire  parameter  set  using  equations  (2.8)-(2.11)  with  the  Xf, 
and  Yf,  calculated  using  equations  (3.3)-(3.4). 

The  LMI  optimization  over  the  second  parameter  subset  Vs,  U  Vd  is  related  with  the  solution  matrices  Xi 
and  Yi  over  the  first  parameter  subset  Vs^  U  Vd-  Thus,  defining  the  order  of  parameter  subsets  may  affect 
the  designed  LPV  controller  over  the  entire  parameter  set. 

Both  methods  I  and  II  require  appropriate  blending  functions  to  blend  solution  matrices  over  the  param¬ 
eter  subset  Vst^  ■  Note  that  an  LPV  controller  provided  by  the  blending  methods  I  and  II  changes  dependently 
on  which  blending  functions  are  selected.  It  is  unknown  how  this  affects  the  closed-loop  performance  of  the 
designed  LPV  controller. 

For  method  III,  blending  matrix  functions  are  calculated  to  minimize  the  induced  norm  7  over  the 
parameter  subset  Vst^  U  Vd-  The  blended  solution  matrices  Vj  and  1),  are  rewritten  as: 

Xbips^Pd)  =  \[X,,{ps)X,{pd)  +  X,{pd)X,,{ps)]  +  ]^[X,,{ps)X2{pd)  +  X2{pd)Xk,{ps%  (3.16) 

yb{ps,Pd)  =  \[YH{ps)Yi{pd)  +  Y,{pd)Y,,{ps)]  +  \[Y,,{ps)Y^{pd)  +  Y^{pd)Y,,{ps)],  (3.17) 

where  “blending  matrix  functions”  Vjj  (p^),  (p^),  l)>j  (Ps),  and  Yf,^  (ps)  are  differentiable  symmetric  matrix 

functions  bounded  over  the  parameter  subset  Vs^  -  To  present  the  blending  matrix  functions,  basis  functions 
9xj{Ps),  hxjiPs),  9vj{Ps),  and  hy.{ps)  for  Xt,^{ps),  Xb^{ps),  Yb^{ps),  and  Yb.,{ps)  are  introduced,  respectively. 
To  satisfy  the  boundary  condition  of  blended  matrices  Xb{ps,Pd)  and  Yi,(ps,pd),  Xi,^(ps),  Xb,{ps),  Yi,^(ps), 
and  Yb,{ps)  are  defined  as: 


where 


XbiiPs)  —  9xo{Ps)I  Y  ^  )  9xj  ^Ps)Xbi. ,  Xb,{ps)  —  hx,(yPs')I  Y  hxj{ps)Xb,, , 


Vbi  (Ps)  —  9yo  {Ps)I  Y  ^  )  9vj  {Ps)Yb^, ,  Yb,  (ps)  —  hyg  {ps)I  Y  ^  )  hy.  {ps)Yb,. , 

i=i  i=i 


9xo(Ps)  =  1,  hx^(ps)  =  0, 
hxM  =  l,  P.,(P.)  =  0, 
p,„(p,)  =  l,  h,,.(p,)  =  0, 


dhxj  (ps)  _  p, 

dhxj  (ps)  _  Q 
()Psi  “ 
dhyj  (Ps)  _  p 
(^Psi 

dhy.  (Ps)  _  p, 


Ps  G  Vsi  —  V, 
Ps  G  Vs,  —  V. 
Ps  G  Vsi  —  V. 
Ps  G  Vs,  —  V. 


(3.18) 

(3.19) 


(3.20) 


Here,  the  basis  functions  are  differentiable  over  the  parameter  subset  Vs^,  and  Xb-^. ,  Xb,. ,  Yb^. ,  and  Yb,. 
are  unknown  constant  matrices  in  77."^".  The  unknown  constant  matrices  can  be  determined  solving  the 
following  LMI  optimization: 


min  7, 

Xh-,  ,  Xh^  ,  Yh-.  ,  Ifco 
^3  ^3  ^3  ^3 


(3.21) 


7 


+  AX, -YAv^)  -  X,CJ,  7“'Si 


CiiX, 

7-iBf 


-In. 

0 


-B 


<0, 


(3.22) 


Ay,  +  Y,At +  Y{2.i^)-CjC2  Y,B^^  7-if7f 

*=i 

Bl^Y, 


-In 


7“'C^i 


0 


0 

-In 


<  0, 


(3.23) 


where 


X,  J-^In 

7“' 4  n 


>0,  Xb>  0,  n  >  0, 


X, 

n 


r^«8 


*=1 


(3.24) 


<7*0^1  +  hxoX2  +  0.5  9xj[XiX,^,  +  Xb^,Xi]  +  Y^  fixAX2Xb2^  +  X,.^^X2\  ^  , 

i=i  k=l 


(3.25) 


9v.yi  +  hyoY2  +  0.5  <7,,.  [  Y,Y,,^  +  Y,,^  Ei  ]  +  ^  4,  [4-4.,  +  4.,  4-]  }  . 


The  procedure  of  designing  an  LPV  controller  over  the  entire  parameter  set  using  this  blending  method 

III  is: 

1.  Partition  the  entire  parameter  set  into  two  subsets  which  have  the  overlapped  subset. 

2.  Solve  the  LMI  optimization  of  equation  (2.3)  for  Xi{p)  and  Yi{p)  over  each  parameter  subset. 

3.  Define  basis  functions  for  blending  matrix  functions  in  equations  (3.18)  and  (3.19),  which  satisfy  the 
boundary  conditions  of  equation  (3.20)  over  each  parameter  subset. 

4.  Solve  the  LMI  optimization  of  equation  (3.21)  over  the  parameter  subset  Vs^  U  Vd- 

5.  Construct  solution  matrices  X  and  Y  over  the  entire  parameter  set  using  the  solution  matrices  Xi, 
X2,Yi,Y2,X,,andY,. 

6.  Construct  an  LPV  controller  over  the  entire  parameter  set  using  equations  (2.8)-(2.11)  based  on 
solution  matrices  X  and  Y  generated  in  step  5. 

In  this  section,  there  were  described  three  methods  to  calculate  the  blended  solution  matrices  X,{ps,pd) 
and  Yb{ps,pd)  which  are  feasible  solutions  of  LMIs  of  equations  (2.4)-(2.5)  over  the  parameter  subset  Vs^  ^Vd- 
Since  all  three  methods  keep  Xi  and  Yi  over  the  parameter  subset  Vs^  —  Vs^ ,  the  blended  LPV  controllers 
generated  by  the  three  methods  are  equal  to  the  regionally  designed  LPV  controller  Ki  over  the  parameter 
subset  Vsi  —  4n  •  Both  methods  I  and  III  can,  also,  keep  the  regional  LPV  controller  K2  over  the  parameter 
subset  'Ps2  —  4n  • 

The  disadvantage  of  method  I  is  that  choosing  blending  functions  is  in  an  ad-hoc  manner.  However, 
method  I  is  the  fastest  to  synthesize  an  LPV  controller  over  the  entire  parameter  set  among  the  three  blending 
methods.  In  method  II,  the  blending  solution  matrices  X,  and  Y,  are  guaranteed  to  be  feasible  solutions 
over  the  parameter  subset  Vs^  ■  This  method  II  requires  solving  the  LMI  optimization  with  the  extra  LMI 
constraints.  In  method  III,  optimal  blending  matrix  functions  are  calculated  to  minimize  the  induced  £.2 
norm  7  over  the  parameter  set  Vs^ ,  based  on  basis  functions  defined  by  a  controller  designer.  It  is  unknown 
how  the  basis  functions  affect  the  blending  matrix  functions  and  the  LPV  control  design. 
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4.  Quasi-LPV  Model  of  F-16  Longitudinal  Axes.  In  this  section,  a  quasi-LPV  model  of  F-16 
longitudinal  axes  is  presented.  The  full  nonlinear  equations  of  an  F-16A  aircraft  are  taken  from  Ref.  [1]. 
The  nonlinear  equations  of  F-16  longitudinal  axes  [1]  are 


qS  cos  a,  ^  ,  c  ^  ^  ^  qS  sin  a,  ^  ,  c  ^  ^  ^ 

V  =  — — — [Cx{a,Se)  +  ■^Cx^{a)q]  -\ - — — [Cz{a,Se)  +  -^Czg{a)q] 


—  g  sin(0  —  a)  +T  cos  a, 
gcos{9  —  a)  sin  a 


qSc 


a  =  -  ^  +  [1  +  ^  (a.  (a)  cos  a  -  (a)  sin  a)l, 

— ^ 

+  —  [Cz(ct,  Se)  COSO!  —  Cx(ct,  Se)  sin a] , 

a  =  ^i^C^Aa)+ACzAa)h  +  f^iC^AaM  + 

V  J-W  C 


e  =  q. 


(4.1) 


(4.2) 

(4.3) 

(4.4) 


Velocity  (V,  ft/sec),  angle  of  attack  {a,  rad),  pitch  rate  (g,  rad/sec),  and  pitch  angle  {9,  rad)  are  the  state 
variables  and  thrust  (T,  lb)  and  elevator  deflection  {Sg,  rad)  are  the  control  variables.  The  aerodynamic 
coefficients  are  lookup  tables  based  on  wind-tunnel  data  from  NASA-Langley  wind-tunnel  tests  on  an  F-16 
aircraft  scaled  model  [16].  The  aerodynamic  data  were  valid  for  a  speed  range  of  100  <  V  <  900  ft/sec, 
an  angle  of  attack  range  of  -10°  <  a  <  45°,  and  an  altitude  range  of  5000  <  h  <  40000  ft.  These  three 
parameters  are  scheduling  parameters  in  the  quasi-LPV  model  derived  for  the  F-16  longitudinal  axes.  Note 
that  V  and  a  are  both  scheduling  parameters  and  states  and  h  is  a  scheduling  parameter  which  enters 
implicitly  into  the  nonlinear  dynamics. 

Unfortunately,  the  control  variable  Sg  does  not  enter  affinely  in  equations  (4.1)-(4.3).  To  derive  a  quasi- 
LPV  model  of  F-16  longitudinal  axes,  it  is  necessary  that  all  controls  be  in  affine  form.  This  is  achieved 
by  transforming  {Sg,  T)  into  synthetic  inputs  Ug  =  [ui  ■  For  details  of  the  transformation,  readers  are 
referred  to  Ref.  [14].  For  the  F-16  quasi-LPV  model,  cos0  and  sin0  are  linearized  about  a  trim  value  9o- 
After  tedious  algebraic  manipulations,  equations  (4.1)-(4.4)  are  rewritten  as 

X  =  A{V,  a,  h)x  +  M{V,  a,  h)us  +  f{V,  a,  h),  (4.5) 

X  =  [V  a  q  9  —  9o]^ , 

where  the  elements  of  matrices  A  €  M  €  ,  and  /  €  are  written  in  Appendix  A. 

Using  the  function  substitution  method  [18,  19],  the  nonlinear  function  f{V,a,h)  can  be  decomposed 
into  quasi-linear  functions  G{V,a,  h)[V  —  Vg  a  —  ag]^  where  G  is  in  Thus,  a  quasi-LPV  model  of  F-16 

aircraft  longitudinal  axes  is  provided.  The  details  of  function  substitution  are  written  in  Appendix  B.  To 
compare  simulation  results  of  the  nonlinear  and  the  quasi-LPV  model  of  the  F-16  aircraft  dynamics,  several 
time  sets  of  inputs  T  and  Sg  are  pre-defined.  In  this  paper,  one  example  of  time  simulations  is  presented  for 
space  limitation.  For  example,  inputs  are  set  as 


=<5 


(4.6) 


To  lb,  0  <  t  <  1,  11  <  t 

T={  -  -  >  - 

Tg  -  2000  lb,  1  <  t  <  11  sec. 


where  Sgg  and  Tg  is  a  trim  value.  The  time  simulation  results  in  Figure  4.1  show  that  the  time  responses 
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Fig.  4.1.  Nonlinear  vs.  quasi-LPV  model  simulations 

of  the  quasi-LPV  and  nonlinear  models  are  almost  matched  to  each  other.  Note  that  the  quasi-LPV  model 
provided  by  the  function  substitution  method  may  change  depending  on  which  one  trim  point  is  selected.  It 
is  unknown  how  this  affects  the  quasi-LPV  model  or  the  LPV  control  design. 

5.  LPV  Control  Problem  Formulation. 

5.1.  Control  Design  Objectives  and  Weighting  Functions.  A  formulation  of  the  LPV  control 
synthesis  of  the  F-16  longitudinal  axes  is  presented  in  this  section.  The  primary  control  objective  for  the 
F-16  longitudinal  flight  controller  is  to  track  velocity  and  flight  path  angle  commands  within  1  ft/sec  and 
0.6°  error  range  in  steady-state  conditions. 

Velocity  and  flight  path  angle  tracking  problems  can  be  formulated  as  model  matching  problems  in  the 
LPV  control  synthesis.  In  this  paper,  we  consider  the  F-16  aircraft  as  an  unmanned  aircraft.  The  ideal 


10 


transfer  function  from  the  flight  path  angle  command  to  the  flight  path  angle  measurement  is  modeled  as 

a  second  order  system,  0 -426(5  +  1.5)  ^  g  rad/sec  bandwidth  and  a  right  hand  zero  at  -1.5  rad/sec. 

For  the  velocity  tracking  problem,  the  ideal  transfer  function  from  the  velocity  command  to  the  velocity 

measurement  is  modeled  as  the  second  order  system  with  0.4  rad/sec  natural  frequency 

-I- 0.8s  -1-0.16  '  ^ 

and  critical  damping. 


Fig.  5.1.  Interconnection  structure  for  the  model  matching  problems. 


A  block  diagram  of  the  interconnection  structure  for  synthesizing  an  LPV  controller  for  the  F-16  longi¬ 
tudinal  axes  is  shown  in  Figure  5.1.  The  velocity,  angle  rate  and  angle  sensors  are  modeled  as  the  first-order 
low  pass  Alters  ^  and  ^  (see  Ref.  [1,  16]).  In  the  LPV  controller  synthesis  models,  the 

sensor  models  shown  in  Figure  5.1  are  approximated  as: 

Sen  =  diag([l,  1,  ^^  ]),  (5.1) 

s  -h  iU 

since  the  interesting  frequency  range  for  LPV  control  synthesis  is  less  than  10  rad/sec.  The  reduced  order 

sensor  models  help  to  reduce  the  overall  state  order  of  an  LPV  controller,  since  the  state  order  of  an  LPV 

controller  is  same  as  that  of  the  augmented  open-loop  system.  The  elevator  actuator  is  modeled  as  the 
90 

first-order  lag  Alter  ^  ^20  limit  is  defined  as  ±  60  deg/sec  in  Ref.  [1,  16].  Since  the  F-16  is  a 

fighter  aircraft,  we  estimated  the  engine  model  as  the  first-order  lag  Alter  ^  ^  ^  which  allows  fast  responses 
in  engine  dynamics.  The  thrust  rate  limit  is  taken  to  be  10000  Ib/sec.  In  the  block  diagram  in  Figure  5.1, 
the  actuator  rate  and  actuator  models  are: 


Act  =  diag([ 


20s  20  4s  4 

s-|-20’s-|-20’s-|-4’s-|-4 


])• 


(5.2) 
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The  performance  weighting  functions  are  chosen  based  on  the  desired  performance  objectives.  The  per¬ 
formance  weighting  function  of  flight  path  angle  W„  ,  ,  is  derived  based  on  the  performance 

^  (s/0.6  +  1) 

objective  to  keep  a  7  tracking  error  less  than  0.6°  for  1  radian  command  in  steady  state  flight.  Since  the 
bandwidth  of  the  ideal  model  from  the  flight  path  angle  command  to  the  flight  path  angle  measurement  is 
0.8  rad/sec,  the  roll  off  frequency  of  the  weighting  function  is  chosen  as  0.6  rad/sec  to  specify  the  tracking 
error  less  than  0.6°  at  the  low  frequency  region  (<  0.6  rad/sec).  The  performance  weighting  function  for  the 
velocity,  ^ ,  is  derived  to  track  velocity  commands  within  1  ft/ sec  error  range  in  steady  state 

flight.  The  unmodeled  dynamics  are  included  in  the  multiplicative  uncertainty  models,  Wmui  and  Wmu2- 
The  uncertainty  weighting  functions  are  rolled  up  in  the  mid-frequency  range  in  order  to  limit  the  bandwidth 
of  the  LPV  controllers.  The  multiplicative  uncertainty  weight  functions  are  set  as 


0.01(s/0.35+  1) 
s/80+1 
0.01(s/0.2+  1) 
s/50  +  1 


The  sensor  noise  models  are  taken  as  constant  across  frequency  to  reduce  the  state  order  of  the  LPV 
controllers.  The  velocity,  angle,  and  angle  rate  sensor  noises  are  modeled  as  white  noises  with  amplitudes  of 
0.8  ft/sec,  0.1°,  and  0.6  deg/sec,  respectively. 

To  solve  LMI  equations  (2.4)-(2.6),  the  basis  function  sets  need  to  be  defined  for  X  and  Y .  There  is 
no  analytical  method  to  choose  the  best  basis  function  set.  Most  often  the  basis  functions  used  are  power 
series[7],  Legendre  polynomials[10],  or  affine  functions  of  scheduling  parameters  [9].  Here,  the  basis  function 
set  for  X{p)  and  Y{p)  is  chosen  as  the  first  order  power  series  {1,  p}  of  the  scheduling  parameter  of  velocity 
to  reduce  the  computation  time  in  the  LPV  control  synthesis.  Note  that  the  basis  functions  of  X  and  Y  do 
not  have  to  be  same. 


5.2.  Blending  Two  Controllers.  In  this  section,  synthesizing  an  LPV  controller  for  the  F-16  longi¬ 
tudinal  axes  using  the  blending  approach  is  demonstrated.  To  apply  the  blending  approach  for  control  of 
the  F-16,  the  entire  parameter  set  (the  flight  envelope)  is  partitioned  into  two  subsets:  high  and  low  altitude 
regions.  Parameter  subsets  are: 


Vi  =  {(V,q;,/i)|100  <  V  <  900  ft/sec, 
V2  =  {(V,q;,/i)|100  <  V  <  900  ft/sec. 


-10  <  a  <45°,  5000  < /i  <  30000  /t}, 
-10  <  a  <  45°,  10000  <h<  40000  ft}. 


Vr\  =  Vi  f\V2- 


(5.3) 


To  use  methods  I  and  II,  blending  functions  are  required  over  the  parameter  subset  Vn,  which  are  satisfied 
with  the  boundary  conditions  of  equations  (3.5)-(3.7).  Blending  functions  hi{h)  and  62 (/i)  are  chosen  as: 


and 


0, 


h  <  10000  ft, 

10000  <h<  30000  ft, 
30000  ft  <  h. 


(5.4) 


bi{h)+b2{h)  =  1. 


(5.5) 
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Using  the  blending  functions,  the  solution  matrices  X,  and  Yi  are  blended  across  the  parameter  subset  Vn 


as: 


Xb(V,a,h)  =  h{h)Xi{V,a)  +  b2{h)X2{V,a),  (5.6) 

Yb{V,a,h)=hi{h)Yi{V,a)+h2{h)Y2{V,a).  (5.7) 


Using  method  I,  an  LPV  controller  Kj  is  constructed  with  the  solution  matrices  X^  and  Y^  using  equa¬ 
tions  (2.8)-(2.11)  over  the  entire  parameter  set  Vi  WP2- 

Using  the  blending  functions  in  equations  (5.4)-(5.5)  and  the  solution  matrices  Xi  and  Yi  over  the 
parameter  subset  Vi ,  the  solution  matrices  X2  and  Y2  are  determined  solving  the  LMI  optimization  over  the 
parameter  set  P2  with  the  extra  LMIs  of  equations  (3.14)-(3.15)  described  in  method  II.  An  LPV  controller 
Kjj  is  constructed  based  on  the  calculated  blended  matrices  Xj,  and  Yj,  using  equations  (2.8)-(2.11)  over  the 
entire  parameter  set. 

To  blend  solution  matrices  V,(U, a)  and  Yi{V,a)  calculated  over  each  parameter  subset,  the  basis  func¬ 
tions  for  blending  matrix  functions  of  equations  (3.18)  and  (3.19)  are  required  to  use  method  III.  In  this 
paper,  the  basis  function  sets  are  chosen  as  {<70,  9i\  for  Vjj  and  ,  and  as  {/iq,  hi}  for  Xb^  and 

respectively. 


9o{h) 


' 

1, 

<  0.5[l  +  cos(i^*7r)], 

0, 

V 


h  <  10000  ft, 

10000  <h<  30000  ft, 
30000  ft  <  h. 


(5.8) 


9i{h) 


0, 

<  0.15[l  +  cos(/^*7r)], 

0, 

V 


h  <  10000  ft, 

10000  <h<  30000  ft, 
30000  ft  <  h. 


(5.9) 


ho{h)  =  1  -  9o{h),  hi{h)  =  -9i{h).  (5.10) 

Note  that  it  is  not  necessary  to  choose  the  same  basis  functions  for  Xb^  and  Yb^ .  From  the  solution  of  the  LMI 
optimization  of  equation  (3.21),  the  matrices  Xb{V,  a,  h)  and  Yb{V,  a,  h)  are  calculated  using  equations  (3.16)- 
(3.19).  An  LPV  controller  Km  is  constructed  from  the  matrices  Xb(V,a,h)  and  Yb{V,a,h)  using  equations 
(2.8)-(2.11)  over  the  entire  parameter  set. 

For  comparison,  an  LPV  controller  Ktot  is  constructed  with  solution  matrices  X  and  Y  over  the  entire 
parameter  set  using  the  conventional  LPV  controller  synthesis  approach.  It  takes  approximately  43  hours  on 
933  MHz  PHI  machine  running  Linux.  Using  the  blending  methods  I,  II,  and  III,  it  takes  approximately  22, 
25,  and  30  hours  on  the  same  machine,  respectively.  The  computation  time  to  synthesize  the  LPV  controller 
for  the  F-16  longitudinal  axes  is  reduced  using  the  blending  approaches. 

6.  Nonlinear  Simnlations.  The  12  state  nonlinear  F-16  aircraft  dynamics  [1]  with  the  synthesized 
LPV  controllers,  Ki,  K2,  Kfot,  Kj,  Kjj,  and  Kjjj,  are  simulated  to  compare  their  time  simulations  in  this 
section.  Recall  that  the  LPV  controllers  Ki  and  K2  are  constructed  over  parameter  subsets  Vi  and  P2, 
respectively.  The  LPV  controller  Kfot  is  synthesized  over  the  entire  parameter  set,  using  the  conventional 
LPV  control  synthesis  methodology.  The  LPV  controllers  Kj,  Kjj,  and  Kjjj  are  the  blended  controller 
using  the  blending  methods  I,  II,  and  III,  respectively.  In  simulations,  the  LPV  controller  is  implemented 
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Fig.  6.1.  Time  simulations  with  the  LPV  eontrollers  around  8000  ft  altitude. 

using  linear  interpolation  at  the  current  values  of  the  scheduling  parameters  between  the  grid  point  solutions. 
The  full  state  models  of  actuators  and  sensors  are  included  in  the  nonlinear  F-16  aircraft  simulations. 

Using  these  LPV  controllers,  the  velocity  and  7  step  responses  are  simulated  around  8000  ft  and  32000 
ft  altitude,  respectively.  The  step  input  sizes  are  a  10  ft/sec  velocity  command  and  a  5°  7  command  at  1  sec. 
The  simulation  results  in  Figure  6.1  show  that  all  velocity  and  7  measurements  match  the  ideal  responses 
within  ±0.25  ft /sec  and  0.06°  tracking  error.  It  is  observed  that  these  LPV  controllers  achieve  the  desired 
performance  objectives. 

The  LPV  controllers  Ki,  Kj,  Kjj,  and  Kjjj  are  exactly  equal  to  each  other  over  the  parameter  subset 
Vi  —Vn,  since  the  blending  methods  I,  II,  and  III  keep  the  regional  LPV  controller  Ki  over  the  parameter 
subset  Vi  —  Vn-  The  step  responses  of  velocity  and  7  with  the  LPV  controller  Kjjj  in  Figure  6.1  represent 
the  step  responses  with  the  LPV  controllers  Ki,  Kj,  and  Kjj  in  the  simulations  around  8000  ft  altitude. 
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The  actuator  deflections  and  their  rates  of  the  time  simulations  are  shown  in  Figure  6.1.  It  is  noted  that  the 
blended  LPV  controller  Km  uses  smaller  actuator  deflections  and  their  rates  than  the  LPV  controller  Ktot 
does  to  achieve  the  performance  objectives. 


h  =  32000  ft 
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Fig.  6.2.  Time  simulations  with  the  LPV  eontrollers  around  32000  ft  altitude. 


The  LPV  controllers  K2,  Kj,  and  Km  are  exactly  equal  to  each  other  over  the  parameter  subset  V2—Vr\- 
The  blended  LPV  controller  Kn  is  different  from  K2  over  the  parameter  subset  1^2— Vn,  since  in  the  method 
II,  the  blended  controller  is  constructed  solving  the  LMI  optimization  of  equations  (2.3)-(2.6)  with  the  extra 
LMI  constraints  of  equations  (3.14)-(3.15).  The  step  responses  with  the  LPV  controllers  Kfot,  Kjjj,  and 
Kjj  are  shown  in  Figure  6.2.  It  is  observed  that  these  LPV  controllers  achieve  the  desired  performance 
objectives  over  the  parameter  set  P2  —Vn-  The  simulation  results  show  that  the  LPV  controller  Kfot  uses 
the  largest  actuator  deflections  and  their  rates  to  achieve  the  performance  objectives. 

For  comparison,  there  is  simulated  one  of  candidate  maneuvers  that  the  F-16  aircraft  flies  from  15000  ft 
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to  32000  ft  across  the  parameter  subset  Vn,  using  20  and  10  ft/sec  velocity  step  and  10°  7  step  commands. 
Velocity,  flight  path  angle,  altitude,  and  angle  of  attack  time  responses  with  the  LPV  controllers  Ktot^  Kjjj, 
Kjj,  and  Kj  are  shown  in  Figure  6.3.  All  the  LPV  controllers  can  achieve  the  performance  objectives  across 
the  parameter  subset  Vn-  It  is  noticed  that  the  velocity  tracking  performance  with  the  blended  controllers 
is  slightly  better  than  the  LPV  controller  Kfot  designed  over  the  entire  parameter  set  Vi  U  P2  using  the 
conventional  LPV  controller  synthesis  methodology.  The  actuator  deflections  and  their  rates  with  the  LPV 
controllers  are  shown  in  Figure  6.4.  The  simulation  results  show  that  the  blended  controllers  use  much 
smaller  actuator  deflections  to  track  the  velocity  and  7  commands. 


Fig.  6.3.  Time  simulations  with  the  LPV  eontrollers  for  the  eandidate  maneuver. 

Sensor  noises  are  integrated  into  the  F-16  aircraft  simulations  for  the  same  situation  of  flying  the  F-16 
aircraft  from  15000  ft  to  32000  ft.  Sensor  noises  for  pitch  rate,  velocity  and  pitch  angle  are  set  as  white 
noises  with  ±0.5  deg/sec,  ±  0.8  ft/sec,  and  ±0.1°  amplitudes.  The  simulation  results  are  omitted  for  space 
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Fig.  6.4.  Actuator  time  responses  with  the  LPV  controllers  for  the  candidate  maneuver. 

limitations.  The  LPV  controllers  Kfot,  Kjjj,  Kjj,  and  Kj  stabilize  the  F-16  aircraft  with  the  sensor  noises 
and  achieve  tracking  performance  objectives  over  the  scheduling  parameter  variations.  The  thrust  signal  and 
its  rate  do  not  exceed  their  limits:  T  <  19000  lb  and  |^|  <  10000  Ib/sec.  Also,  the  elevator  actuator  and 
its  rates  do  not  exceed  their  limits:  |(5e|  <  25°  and  |^|  <60  deg/sec.  The  blended  controllers  Km  and  Kj 
can  preserve  the  performance  level  of  the  regional  LPV  controllers  Ki  and  over  each  parameter  subset 
Vi  and  V-i  ■ 

Note  that  method  I  requires  the  slow  parameter-varying  blending  functions.  In  this  example,  the  param¬ 
eter  intersection  range  is  wide  enough  to  provide  the  slow-varying  blending  functions  for  method  1.  When 
the  parameter  intersection  space  is  narrow,  method  III  is  appropriate  to  blend  two  LPV  controllers  over  the 
parameter  intersection  space. 
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7.  Summary  and  Discussion.  In  this  paper,  the  methodologies  of  blending  LPV  controllers  are 
discussed  to  preserve  the  performance  level  over  each  parameter  subset  and  to  reduce  computation  time  to 
synthesize  an  LPV  controller  over  the  entire  parameter  set.  The  blending  approaches  are  to  design  each 
LPV  controller  for  the  set  of  a  small  number  of  scheduling  parameters  over  the  parameter  subsets  and  blend 
all  controllers  scheduled  on  the  entire  scheduling  parameters  using  the  blending  functions. 

The  quasi-LPV  model  of  the  F-16  longitudinal  axes  is  provided  by  a  function  substitution  method  over 
the  entire  flight  envelope  to  facilitate  synthesis  of  an  LPV  controller.  The  two  LPV  controllers  of  the  F-16 
longitudinal  axes  are  synthesized  as  functions  of  velocity  and  angle  of  attack  at  two  regions:  low  and  high 
altitudes,  respectively.  The  two  LPV  controllers  are  blended  into  a  single  LPV  controller  as  functions  of 
velocity,  angle  of  attack,  and  altitude  over  the  entire  flight  envelope,  using  the  three  blending  approaches.  It 
is  noted  from  the  nonlinear  time  simulations  that  the  blended  controllers  achieve  better  performance  than 
the  LPV  controller  constructed  using  the  conventional  LPV  synthesis  methodology.  The  blended  controllers 
constructed  using  methods  I,  II,  and  III  achieve  the  performance  objectives  and  stabilize  the  closed-loop 
system  of  the  F-16  aircraft  with  sensor  noises. 
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Appendix  A.  Nonlinear  Eqnations  of  F-16  Longitndinal  Axes.  Prom  nonlinear  equations  (4.1)- 
(4.4),  the  nonlinear  equations  of  the  P-16  longitudinal  axes  are  rewritten  as  follows: 
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The  gain  matrix  M  is  constructed  as  lookup  tables  and  Ug  are  synthetic  inputs.  The  term  of  M{V,a,h)u, 
can  represent  the  following  terms: 


M{V,a,h)ug{V,a,Se,T) 


f  [Cz, sin  a +  Cx{a,Sg)  cos  a] +  ^-^T 
■:^[CzJe  cosQ:-C'x(Q;,(5e)sinQ:]  - 

^^CM(a,S,)  +  ^ACzJe 

^yy  ^yy 

0 


The  detailed  methods  to  determine  M  and  synthetic  inputs  Ug  are  available  in  Ref.  [14].  The  synthetic 
inputs  can  vary  in  the  range  of  — 1  <  ui  <  1  and  0  <  U2  <  1.  The  units  of  ui  and  U2  are  25°  and  19000  Ih, 
respectively. 

Appendix  B.  Decomposition.  Set  the  state  variables  in  a  quasi-LPV  model  of  the  F-16  aircraft 
dynamics  as  R  —  Vg,  a  —  ag,  q  —  Qg  and  9  —  0g,  where  Vg,  ag,  qg  and  9g  represent  a  trim  point.  Then, 
equation  (A.l)  is  rewritten  as 
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and 


^1=2:  -2:0,  ^2=W-Wg,  U  =  Ug-Ug, 


z  =  [V  a]’^ ,  w  =  [q  9]^ . 
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To  provide  a  quasi-LPV  model  of  F-16  aircraft  dynamics  from  equation  (B.l),  the  term  F{z)  should  be 
decomposed  into  linear  parameter  varying  functions  written  as 


F(2)=F(2o+6) 
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where  gz  G  and  G  .  There  are  an  infinite  number  of  possible  solutions  of  gz  and  g^  to  satisfy 

equation  (B.4).  To  determine  functions  gz  and  g^,  more  constraints  are  required.  In  this  paper,  the  variation 
of  gz  and  g^  over  the  entire  flight  envelope  is  minimized.  With  these  constraints,  an  optimization  problem  is 
formulated  to  determine  gz  and  g^.  For  example,  to  determine  the  first  row  of  gz,  an  optimization  problem 
is  follows: 
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where  Fi  is  the  first  row  of  F.  To  make  the  optimization  problem  of  equation  (B.5)  computationally  tractable, 
the  continuous  constraints  are  evaluated  at  grid  points  over  the  parameter  set.  Thus,  the  matrix  can  be 
determined  at  every  grid  point  of  Using  solutions  of  the  optimization  for  g^  and  g^,  a  quasi-LPV  model 
of  F-16  longitudinal  axes  is  written  as: 
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